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Abstract 

We study the critical properties of three-dimensional O(iV) models, for N = 2, 3, 4. 
Parameterizing the leading corrections-to-scaling for the i] exponent, we obtain a re- 
liable infinite volume extrapolation, incompatible with previous Monte Carlo values, 
but in agreement with e-expansions. We also measure the critical exponent related 
with the tensorial magnetization as well as the v exponents and critical couplings. 
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1 Introduction 



The study of continuous spin models in low dimensions has been very useful 
either for the knowledge of some physical systems directly associated (mainly 
in condensed matter), or as toy systems for studying relativistic field theories. 

The 0(4) model in three dimensions has been conjectured to be on the same 
universality class as the finite-temperature chiral phase transition of QCD 
with massless flavors [1]. The 0(4) universality class has also appeared in 
perturbation-theory studies of spin-systems for which the 0(3) symmetry of 
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the action is fully broken on the low temperature phase [2]. It is well known 
that the 0(3) model in two dimensions offers a play-ground to explore asymp- 
totic freedom [3], although there has been recently some controversy on this 
point [4]. Regarding the applications to condensed matter physics, let us re- 
mind that the three-dimensional 0(3) model is the low-temperature's effective- 
model for a bidimensional quantum antiferromagnet [5]. It also appears as a 
limiting case of a Z2-gauge lattice model for nematic phase transitions of liquid 
crystals [6]. Finally, the 0(2) model in three dimensions is known to be in the 
same universality class as superfluid 4 He. 

The 0(N) nonlinear sigma models in three dimensions have been extensively 
studied either with analytical [7,8] or numerical methods, obtaining very ac- 
curate results in the determination of the critical properties. In particular the 
critical exponents have been measured for N < 4 with a precision greater than 
a 1% [9-11]. However, the finite-size effects have not been considered in a sys- 
tematic way. We find this procedure not harmful for v exponent measures, but 
quite dangerous for r\ exponents determinations. 

Another point we focus on, is the study of the second-rank tensorial mag- 
netization. This composite operator is needed as an order-parameter in some 
applications, for instance when the usual magnetization vanishes, as in studies 
of nematics [6]. It is also the relevant order parameter for systems where the 
0(3) symmetry is fully broken [2]. Naively one could expect that the tenso- 
rial magnetization scales just as the square of the vectorial one, that is, /?t is 
twice (3, as a mean field calculation predicts (it was also assumed in ref. [6]). 
Nevertheless we shall show that this is not so. A somehow related question, 
which has aroused interest lately, is the presence of some pathologies in the 
renormalization of composite operators in 2 + e expansions [12]. 

In this paper we consider the 0(2), 0(3) and 0(4) models, centering mainly in 
the measures of magnetic exponents using a finite-size scaling method which 
is specially useful for observables that change rapidly at the critical point. 
In addition to the standard magnetic exponents we measure those related 
with tensorial excitations. The critical couplings and v exponents are also 
studied. We are specially interested in the measure of v in the 0(4) case 
in order to compare with the results obtained in the antiferromagnetic RP 2 
model [13], which is related with 0(4) when there is a total breakdown of its 
0(3) symmetry. 
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2 The model 



We consider the usual Hamiltonian 

n = -p ]T Vi ■ v v (i) 

<i,3> 

where V{ is a N components normalized vector, and the sum is extended over 
first neighbor pairs. 

It is well known that this model undergoes a second order phase transition 
for which the normalized magnetization M = yY,i v i is an order parame- 
ter (V is the lattice volume). As this model can be simulated using cluster 
algorithms [14], it is possible to thermalize very large lattices. 

We are also interested in studying the behavior of composite operators that 
could be related with bound states from a quantum field theory point of view. 
One can construct orthogonal states to that generated by the fundamental 
field, just ensuring that the composite operator transforms as a higher order 
irreducible representation. The simplest representation beyond the fundamen- 
tal one is the second rank tensorial representation. The associated tensorial 
magnetization can be written as 

A*°* = ^E(^?-^*)- ( 2 ) 

i 

As it happens with the vectorial magnetization, the mean value of M. is zero 
in a finite lattice, so, in the Monte Carlo simulation we have to construct 
an estimator that avoids the tunneling effects. We define the (normalized) 
magnetizations as 




We also define the associated susceptibilities as 

X = V(M 2 ) , Xr = v(trM 2 ). (4) 

The critical behavior of those quantities is expected to be 

Moct^ , M T cxt pT , xoc|t|~ 7 , XT0c|t|~ 7T , (5) 

where t is the reduced temperature. We expect the tensorial exponents to be 
also related through the (hyper) scaling relation 7t + 2/?t = vd. 
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3 The method 



To compute the critical exponents we have used a finite-size scaling analysis. 
Specifically, we have used the method of refs. [13] that consists in the com- 
parison of observables on several pairs of lattices at the same coupling. The 
mean value for the operator O, measured in a length L lattice, at a coupling 
(3 in the critical region, is expected to behave as 

(0(L, (3)) = L*°l" \Fo(Z(L, (3)/L) + L-"G Q m, P)/L) + ...], (6) 



where Fo and Go are smooth scaling functions for this operator, £ is a measure 
of the correlation length and xo is the critical exponent associated with O. 
Let us recall that uj is an universal exponent related with the first irrelevant 
operator. The dots stand for further scaling corrections. We have also dropped 
a £ _a; term, negligible in the critical region. To eliminate the scaling functions 
we compare the measures in two different lattice sizes (Li,L 2 ) at the same 
coupling. Let be 

_ (Q(L 2 ,/3)) _ /v F {^L 2 ^)/L 2 ) 
Qo ~WMM- S F muP)/L 1 ) +0{L } ' (7) 



where s — L 2 /L ± . It is easy to find the coupling value where = s. Measuring 
O at that point we obtain the critical exponent from 

Qo\ Qi=s = s*°/» + 0(L-"). (8) 



We remark that even if O is a fast varying function of the coupling in the 
critical region, as the magnetization, the statistical correlation between Qo 
and allows a very precise determination of the exponent. 

For the correlation length, we use a second momentum definition [15] which 
is easy to measure and permits to obtain an accurate value: 

( ( x/F-i y /2 

4 U S in 2 (7r/iW ' W 



where F is defined as the Fourier transform of the two-point correlation func- 
tion at minimal momentum ((2tt/L, 0, 0) and permutations). 
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4 Critical exponents 



For the Monte Carlo simulation, we have used the Wolff's embedding algo- 
rithm with a single cluster update [14]. We have simulated in lattice sizes 
from L = 8 to L = 64 on the critical coupling reported on refs. [9-11]. We 
have used the spectral density method to extrapolate to the neighborhood of 
these couplings. We have updated 25 million clusters for 0(2) and 0(3) and 
50 million in the 0(4) case. The autocorrelation times are very small in all 
cases (not larger than a hundred of clusters). The runs have been performed 
on several workstations. 

We have used the operator d£/ d/3 to obtain the v exponent (x d ^ d/3 = 1 + v). 
For the magnetic exponents, we use the total magnetization (x M = ($) as 
well as the corresponding susceptibility (x x = 7). From them we obtain the 
77 exponent using the scaling relations 7/V = 2 — r),2(3/v = d — 2 + rj. For 
the tensorial channel we obtain a different set of exponents that we denote as 
/3t, 7t, Vt respectively. 

In table 1 we report the results for the exponents displaying also the used 
operator. We have checked other observables as well as other definitions of 
the correlation length but, in all cases, either the corrections-to-scaling or the 
statistical errors are greater. 



5 Infinite volume extrapolation 

In table 1 corrections-to-scaling are clearly visible for the rj exponent. To con- 
trol them, we need an estimation of the u exponent. To measure this index and 
the critical coupling [13], we study the crossing between the Binder cumulant 
of the magnetization for different lattice sizes, as well as the corresponding for 
£/L. The shift from the critical coupling of the crossing point for lattices of 
sizes Li and L 2 behaves as [16] 

A/3 Ll ' i2 oc 1 ~ S - L^~». (10) 
s» — 1 

As the crossing point for the Binder cumulant and £/L tends to the critical 
coupling from opposite sides, it is convenient to fit all the data together. We 
have carried out two types of fits, one fixing the smaller lattice and the other 
fixing the L 2 /Li ratio. In table 2 we present the results obtained using the 
full covariance matrix, for L\ — 8 and s = 2. 
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Table 1 

Critical exponents obtained from a finite-size scaling analysis using data from lat- 
tices of sizes L and 2L for the O(N) models. In the second row we show the operator 
used for each column. 







V 


V 


VT 


N 


L 


d^/d(5 


X 


M 


XT 


M T 


2 


8 


0.683(3) 


0.0252(10) 


0.0297(11) 


1.499(2) 


1.506(2) 




12 


0.678(3) 


0.0303(11) 


0.0333(13) 


1.496(3) 


1.500(3) 




16 


0.672(3) 


0.0329(12) 


0.0355(12) 


1.494(2) 


1.497(3) 




24 


0.676(4) 


0.0344(12) 


0.0355(13) 


1.494(3) 


1.495(3) 




32 


0.670(3) 


0.0366(12) 


0.0387(13) 


1.494(3) 


1.496(3) 


3 


8 


0.724(3) 


0.0300(10) 


0.0317(10) 


1.432(2) 


1.437(2) 




12 


0.712(4) 


0.0337(9) 


0.0352(10) 


1.4312(17) 


1.4342(17) 




16 


0.712(4) 


0.0344(11) 


0.0354(12) 


1.428(2) 


1.429(2) 




24 


0.716(5) 


0.0378(12) 


0.0385(13) 


1.4320(18) 


1.4335(18) 




32 


0.711(5) 


0.0371(11) 


0.0377(12) 


1.428(2) 


1.430(2) 


4 


8 


0.752(2) 


0.0307(6) 


0.0316(6) 


1.3735(10) 


1.3767(11) 




12 


0.747(3) 


0.0338(5) 


0.0345(6) 


1.3771(10) 


1.3790(10) 




16 


0.754(4) 


0.0341(7) 


0.0344(7) 


1.3770(14) 


1.3781(13) 




24 


0.757(4) 


0.0348(5) 


0.0349(5) 


1.3771(10) 


1.3774(11) 




32 


0.753(5) 


0.0359(9) 


0.0361(10) 


1.3753(18) 


1.3759(18) 



We observe a value for the uj exponent compatible with the 0.78(2) value 
obtained from (^-expansions for the 0(2) and 0(3) models [7]. In the 0(4) 
case the result is almost twice. Lacking a theoretical prediction, this could be 
interpreted in two ways, the exponent could truly be so big, or it might be 
that the coefficient of the leading corrections-to-scaling term is exceedingly 
small. 

The values obtained for the critical couplings are extremely precise, and com- 
patible with the most accurate previous determinations by Monte Carlo sim- 
ulations [17,10,11]. Let us remark that the conjectured value for the critical 
coupling of the 0(3) model [18], (3 C = log 2, is ten standard deviations away 
from our measures. 

To control finite-size effects the most common procedure is to pick the smaller 
lattice for which there is no scaling corrections. That is, one uses a log-log 
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Table 2 

Fits for P c (oo) and the corrections-to-scaling exponent u>. The second error bar in 
co is due to the variation of v exponent within a 1% interval (the size of previously- 
published error bars). 



N 


Fit 


xVd-o.f. 


UJ 


&(°°) 


2 


Li = 8 


11 /8 


0.86(12)(3) 


0.454169(4) 


s = 2 


9.8/6 


0.81(12)(1) 


0.454165(4) 


3 


Li = 8 


2.0/8 


0.64(13)(2) 


0.693001(10) 


s = 2 


2.3/6 


0.71(15)(1) 


0.693002(12) 


4 


Li = 8 


11.9/8 


1.80(18)(6) 


0.935858(8) 


s = 2 


8.0/6 


1.85(21)(2) 


0.935861(8) 



plot, discard small lattices, and stop when the value of critical index stabilizes. 
Using this method, it has been determined that for the 0(2) model it is enough 
to use lattices with sizes L > 16 [9], L > 12 for the 0(3) model [10], and L > 8 
for the 0(4) model [11]. We shall confirm this assumption for the v exponent 
but no for the magnetic ones. 

It might not be possible to find a safe L min lattice (in fact, to find scaling 
corrections is just a matter of statistical accuracy), we thus need an extrap- 
olation procedure. With an u estimation, we can extrapolate to the infinite 
volume limit, with an ansatz for the exponent xq- 



xq 

V 



Xq 
V 



oc L~ u . 



(L,2L) 



The situation is fairly different for the v and r/ exponents, therefore we shall 
discuss them separately. 



5.1 rj type exponents. 



Due to the high accuracy that we get from the statistical correlation between 
the measures of £ and \i we are & ble to resolve finite-size corrections. One 
could wonder if only the first correction term (oc L~ u ) is needed. To check 
it we have used an objective criterium: we perform the fit from L min , then 
we repeat it discarding L min and check if both fitted parameters (slope and 
extrapolation) are compatible. If that is the case, we take the central value 
from the fit with L min and the error bar from the fit without it. In the three 
models, we find that L min = 8 is enough for this purpose. 

For the 0(2) and 0(3) models, we have a very precise knowledge of uj = 
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0.78(2), from series analysis. We plot in figure 1 our data as a function of L w . 
For 0(2), the fit has x 2 /d.o.f. = 0.85/3. For 0(3) we obtain x 2 /d.o.f. = 2.98/3. 
We then find no reason to expect higher order corrections to be significant. 
We have also performed a simulation on a L = 6 lattice, finding that, for the 
0(2) model, the corresponding r\ value is one standard deviation away from 
the corresponding point in the fit performed for L m i n = 12. Therefore, it seems 
even reasonable to keep the error from the L m i n = 8 fit (rj = 0.0424(14)). For 
the 0(3) model, the L = 6 point is three standard deviations away from the 
fit for L min = 12 due to higher order corrections-to-scaling. Notice in fig. 2 
that the slope for 0(2) is significantly larger, which could mask higher order 
corrections. For the 0(4) model, lacking a theoretical knowledge of ou, there 
are stronger uncertainties. We find the change from uo = 0.78(2) for 0(2) and 
0(3), to u = 1.8(2) for 0(4) very surprising (see table 2). This might arise 
from an unexpected cancelation of first order scaling-corrections for the Binder 
cumulant and the correlation length. We show in table 3 the extrapolated 
values for 77 exponents. For the 0(4) model we present the extrapolation with 
a = 0.78 (xVd.o.f. = 2.9/3) and with u = 1.8(2) with x 2 /d.o.f. = 1.0/3. 
Both values are hardly compatible. 

The determination of rj from x an d M are of course coincident in all cases 
due to the strong statistical correlation between both observables. We find 
that the value for 0(2) and 0(3) are compatible with e-expansions and not 
too far from g-expansions, but incompatible with the previous Monte Carlo 
results [9-11]. For 0(4) both rj values are significantly lower than for the other 
models which is not surprising because, in the large N limit, this exponent 
should go to zero. 

Let us consider the t]t exponent. From table 1 corrections-to-scaling are not 
self-evident. If one adopts this optimistic point of view and performs a mean 
of all values in the table one finds that only for L > 16, x 2 /d.o.f. becomes ac- 

Table 3 

Infinite volume extrapolation for i] from \ an< l M. The second error is due to the 
uncertainty on to. For the 0(4) model, the first column has been calculated with co = 
0.78 and the second with u> = 1.8(2). In the last three rows we present, respectively, 
the results for this exponent from e, g-expansions and previously published Monte 
Carlo values. 





0(2) 


0(3) 


0(4) 


X 


0.0424(23)(2) 


0.0413(15)(1) 


0.0384(12) 


0.0359(6)(3) 


M 


0.0421(25)(2) 


0.0414(18)(1) 


0.0381(13) 


0.0359(6)(2) 


e-expansion 


0.040(3) [7] 


0.040(3) [7] 


0.03(1) [20] 


g-expansion 


0.033(4) [8] 


0.033(4) [8] 




Previous MC 


0.024(6) [9] 


0.028(2) [10] 


0.025(4) [11] 



8 



0.045 




fr 0.035 



0.030 



0.025 



Fig. 1. r] estimation from \ f° r pairs of lattices of sizes L and 2L, as a function of 
L~ u , lo = 0.78, in the 0(2) model (upper part), and 0(3) model (lower part). Lines 
correspond to the L m i n = 12 fit. 

ceptable. Errors for the means are much smaller than for individual measures, 
and so, maybe smaller than scaling-corrections. We show these mean values 
on table 4. We think safer to perform a fit to the functional form (11). We 
show the extrapolated values on table 4. Errors coming from the uncertainty 
on u are smaller than a 10% of the quoted errors. Even more, the two values 
for 0(4) are compatible. The value of x 2 /d.o.f. for 0(4) and 0(2) models is 
about 0.6/3. For the 0(3) model is higher (about 6/3). Regarding the conjec- 
tured relation 2(3 = (3t proposed in ref. [6], notice that it can be formulated 
equivalently as 1 + 2i] = t\t which is absolutely ruled out by our data. 

For the comparison between the critical exponents of the 0(4) model and 
those of the antiferromagnetic RP 2 one, notice that the 77 value for 0(4) is fully 
compatible with the corresponding exponent for the staggered magnetization 
of the RP 2 model( 77 stag = 0.0380(26) [13]). The r/ exponent associated with 
the usual magnetization for the RP 2 model (r] RP 2 = 1.339(10) [13]) is not 
compatible with r] T , the difference being of order r\. 



5.2 v exponent. 



In figure 2, we plot the estimation of v that we get from pairs of lattices 
of sizes L and 2L as a function of L~ w . Indeed one would be tempted to 
claim that for the 0(4) model, finite-size effects are beyond our resolution for 
L > 8, as stated in ref. [11]. For the 0(3) model this also seems to be the 
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Fig. 2. v estimation from pairs of lattices of sizes L and 2L, as a function of L 78 , 
in the 0(4) (upper part), 0(3) (middle part), and 0(2) model (lower part). 

case for L > 12 [10]. However the question is not so clear for the 0(2) model. 
Nevertheless, the experimental value is known to be z/ exp = 0.6705(6) [19], and 
we do find that for L > 16 our v estimation is consistent with it. If we actually 
believe that the corrections-to-scaling are negligible for L > L min , we can take 
the mean of the safe lattices, getting 



i/ (2) = 0.6721(13) ,r7ci-cf. = 1.7/2, 
u 0(3) = 0.7128(14) ,x 2 /d.o.f. =0.6/3, 
z/ 0(4) = 0.7525(10) ,x 2 /d.o.f. = 3.4/4. 



(12) 



Error bars decrease strongly compared to the data in table 1, therefore it is 
not clear if scaling-corrections are still negligible. A more conservative point of 
view would ask for the consideration of these corrections. For this we use the 
ansatz (11). However, from the plot in figure 2, it seems clear that the scaling 
corrections for L — 8, 12 in the 0(2) model, and L = 8 in the 0(3) model, 
could hardly be linear on L~ u . Performing the fit suggested by equation (11), 
with u = 0.78(2) for 0(2) and 0(3) [7] from the safe L min , we get a compatible 

Table 4 

For every model we show in the first column the mean values for tjt from \T and 
Mt for L > 16. The second columns display the infinite volume extrapolation for 
rjT from xt and My. For the 0(4) model, the second column has been calculated 
with uj = 0.78 and the third with to = 1.8(2). 





0(2) 


0(3) 


0(4) 


XT 


1.494(1) 


1.489(4) 


1.431(1) 


1.427(3) 


1.3766(5) 


1.374(5) 


1.376(2) 


M T 


1.496(1) 


1.491(4) 


1.429(1) 


1.427(3) 


1.3773(5) 


1.375(5) 


1.376(2) 
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value, but with a very increased error bar: 



^o(2) = 0.670(10), 
z/ 0( 3) = 0.711(10). 

In the 0(4) case, we perform the extrapolation with both uj values. The extrap- 
olation does not increase significantly the error bars. Therefore we consider 
two values of L m i n . For L m i n = 8 we obtain v = 0.754(3) for uj = 0.78 and 
v = 0.7531(15)(3) for uj = 1.8(2). On the other hand, with L min = 12 we get 
v = 0.765(8) for u = 0.78 and v = 0.7585(34)(7), for to = 1.8(2). 

Regarding the relation with the antiferromagnetic RP 2 model {y = 0.783(11)), 
we find it unlikely that both values could coincide, but we cannot rule it out. 
Notice that if we accept the value uj = 1.8(2) we would find a significant 
difference on this universal exponent from RP 2 (uj = 0.85(5)). 



6 Conclusions 



We have obtained accurate measures of critical exponents and couplings for 
three dimensional O(N) models. The method used, based on the finite-size 
scaling ansatz, has a remarkable performance when computing magnetic ex- 
ponents. The values for 77 exponents presented are incompatible with previous 
Monte Carlo results and they have smaller statistical errors. We should point 
out that our values agree with those obtained with 4 — e expansions, and are 
not incompatible with those computed by means of g expansions, in opposition 
with previous results. We also show that the statistical accuracy that can be 
reached on three dimensional systems is such that the strongest uncertainties 
come from finite-size effects, for which a reliable theoretical parameterization 
would be highly desirable. 

We present measures of the exponent corresponding to the tensorial magne- 
tization, which is not twice the corresponding to the usual magnetization as 
previously stated. 

We increase significantly the precision of previous measures of the exponents 
of 0(4), showing that the differences with the v exponent from the values 
of the antiferromagnetic RP 2 model [13] are two standard deviations appart. 
However whether they belong to the same universality class or not is not yet 
completely established. 

This work has been partially supported by CICyT AEN93-0604-C03-02 and 
AEN95-1284-E. 



11 



References 



[I] R. Pisarski and F. Wilczeck, Phys. Rev. D29 (1984) 338; F. Wilczek, 
Int. J. Mod. Phys. A7 (1992) 3911; K. Rajagopal and F. Wilczek, Nucl. Phys. 
B399 (1993) 395. 

[2] P. Azaria, B. Delamotte and T. Jolicoeur, Phys. Rev. Lett. 64 (1990) 3175; P. 
Azaria, B. Delamotte, F. Delduc and T. Jolicoeur, Nucl. Phys. B408 (1993) 
485. 

[3] A. M. Polyakov, Phys. Lett. B59 (1975) 79. 

[4] S. Caracciolo, R. G. Edwards, A. Pelissetto and A. D. Sokal, Phys. Rev. Lett. 
76 (1996) 1179; A. Patrascioiu and E. Seiler, Phys. Rev. Lett 76 (1996) 1178. 

[5] E. Manoussakis and R. Salvador, Phys. Rev. B40 (1989) 2205. 

[6] P. E. Lammert, D. S. Rokhsar and J. Toner, Phys. Rev. Lett. 70 (1993) 1650. 

[7] J. C. Le Guillou and J. Zinn-Justin, Phys. Rev. B21 (1980) 3976. 

[8] J. C. Le Guillou and J. Zinn-Justin, Nucl. Phys B257[FS14] (1985) 867; 
J. C. Le Guillou and J. Zinn-Justin, J. Physique Lett. 46, (1985) L137 

[9] A. P. Gottlob and M. Hasenbusch, Physica A201 (1993) 593. 

[10] C. Holm and W. Janke, Phys. Lett. A173 (1993) 8. 

[II] K. Kanaya and S. Kaya, Phys. Rev. D51 (1995) 24. 
[12] F.J. Wegner, Z. Phys. B78 (1990) 33. 

[13] H. G. Ballesteros, L. A. Fernandez, V. Martin-Mayor, A. Munoz Sudupe. 

Preprint |hep-lat/9511003| . To appear in Phys. Lett. B; Preprint [hep-| 
| lat/9605037[ 

[14] U. Wolff, Phys. Rev. Lett. 62 (1989) 3834. 

[15] F. Cooper, B. Freedman and D. Preston, Nucl. Phys. B210 (1989) 210. 
[16] K. Binder, Z. Phys. B43 (1981) 119. 

[17] A. P. Gottlob and M. Hasenbusch, J. Stat. Phys. 77 (1994) 919. 

[18] R. G. Brown and M. Ciftan, Phys. Rev. Lett. 76 (1996) 1352. 

[19] L. S. Goldner and G. Ahlers, Phys. Rev. B45 (1992) 13129. 

[20] G. Baker, D. Meiron, and B. Nickel, Phys. Rev. B17 (1978) 1365; Compilation of 
2-pt. and J^-pt. graphs for continuous spin models, University of Guelph Report, 
1977 (unpublished). 



12 



